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We have studied the extended Hubbard model with pair hopping in the atomic limit for arbitrary 
electron density and chemical potential. The Hamiltonian considered consists of (i) the effective 
on-site interaction U and (ii) the intersite charge exchange interactions /, determining the hopping 
of electron pairs between nearest-neighbour sites. The model can be treated as a simple effective 
model of a superconductor with very short coherence length in which electrons are localized and only 
electron pairs have possibility of transferring. The phase diagrams and thermodynamic properties of 
this model have been determined within the variational approach, which treats the on-site interaction 
term exactly and the intersite interactions within the mean-field approximation. We have also 
obtained rigorous results for a linear chain (d = 1) in the ground state. Moreover, at T = some 
results derived within the random phase approximation (and the spin-wave approximation) for 
d — 2 and d — 3 lattices and within the low density expansions for d = 3 lattices are presented. Our 
investigation of the general case (as a function of the electron concentration n and as a function of 
the chemical potential fi) shows that, depending on the values of interaction parameters, the system 
can exhibit not only the homogeneous phases: superconducting (SS) and nonordered (NO), but also 
the phase separated states (PS: SS-NO). The system considered exhibits interesting multicritical 
behaviour including tricritical points. 

PACS numbers: 71.10.Fd, 71.10.-w, 74.20.-z, 74.81.-g, 64.75. Gh 
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I. INTRODUCTION 

There has been much interest in superconductivity 
with very short coherence length. This interest is due to 
its possible relevance to high temperature superconduc- 
tors (cuprates, doped bismuthates, iron-based systems, 
fuUerenes) and also to the several other exotic supercon- 
ducting materials (for a review, see [1, 2_ and references 
therein). It can also give relevant insight into behaviour 
of strongly bounded fermion pairs on the optical lattices. 

The phase separations involving superconducting 
states have been evidenced in a broad range of currently 
intensely investigated materials including iron-pnictides, 
cuprates and organic conductors (see for example |2HI2 
and references therein). 

In our work we will study a model which is a simple 
generahzation of the standard model of a local pair su- 
perconductor with on-site pairing (i. e. the model of hard 
core bosons on a lattice [U fT3lll5| ) to the case of finite 
pair binding energy. The Hamiltonian considered has the 
following form: 



H = VY. 



n-i^nii' 



(iJ) 



{ptp~, 



Pj P^ 



/iVn, 



(1) 



C+C+- 



where fii = J2a'^i<^- ^i<y = ^ta^i-^i Pi = (Pi )^ 
Cicr (c^) denotes the annihilation (creation) operator of 
an electron with spin a =t,4- at the site i, which satisfy 
canonical anticommutation relations 



I'^'icr 5 C..(j/ / — ^ij^aa'j \Ciai^j(y'i — 1 



^ia-i ^ja' 



} = 0, (2) 



where 5ij is the Kronecker delta. 'Y^u j) indicates the 
sum over nearest-neighbour sites i and j independently. 
z will denote the number of nearest neighbours, fi is 
the chemical potential, depending on the concentration 
of electrons: 






(3) 



with < n < 2 and N is the total number of lattice 
sites. There are two interaction parameters of the model: 
(i) the intersite pair hopping interaction /, determining 
the electron pair mobility and responsible for the long- 
range superconducting order in the system and (ii) the 
on-site density-density interaction U, which contributes 
(together with /) to the pair binding energy by reducing 
{U > 0) or enhancing ([/ < 0) its value. To simplify our 
analysis we do not include in Hamiltonian (II]) the single 
electron hopping term (t^) as well as other inter-site in- 
teraction terms. This assumption corresponds to the sit- 
uation when the single particle mobility is much smaller 
than the pair mobility and can be neglected. In the pres- 
ence of the tij term, model (ITl is called the Penson-Kolb- 
Hubbard (PKH) model [16-23 and the effects of Uj ^ 
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will be pointed out in section VI 

Model (fTl) can be considered as a relatively simple, 
effective model of a superconductor with local electron 
pairing. Moreover, the knowledge of the tij = limit can 
be used as starting point for a perturbation expansion 
in powers of the hopping tij and provides a benchmark 
for various approximate approaches analyzing the corre- 
sponding finite bandwidth models. 

In terms of pseudospin charge operators {pi} Hamil- 
tonian (IT|) takes the form of a special XY-model with 
a single ion anisotropy in a transverse field (in the z- 



direction) 

H = 2C/^(pD' - / E (Ptpj + PtPi) + 

-A2E(2/5,- + l)-iVf//2, (4) 

i 

where pi = \{ni — 1) and p, — ji — C//2 corresponds to 

the external magnetic field. The operators {pi} satisfy 
standard commutation relations for momentum opera- 
tors: 

[plpf] = ±pfk„ [Ptrp-J] = 2p^<5„-, (5) 

what follows directly from ^. However, one should 
stress that {/5f} takes on four values {—1/2,0, 0, +1/2}, 
as for S = 1 model with doubly degenerated value, 
whereas {p[)'^ = {pT)"^ =0, as in the S — 1/2 system. 
Condition dsl for the electron concentration reads: 



--(n-l) 



^E(^^ 



(6) 



It means that magnetic field p. is determined by a fixed 
value of magnetization (n — l)/2. 

The interactions U and / will be treated as the effec- 
tive ones and will be assumed to include all the possi- 
ble contributions and renormalizations like those com- 
ing from the strong electron-phonon coupling or from 
the coupling between electrons and other electronic sub- 
systems in solid or chemical complexes [T. In such 
a general case arbitrary values and signs of U and / 
are important to consider. It is notable that formally 
/ is one of the off-diagonal terms of the Coulomb in- 
teraction lij — —{\/2){ii\e^/r\jj) [25) describing a part 
of the so-called bond-charge interaction, and it sign of 
the Coulomb-driven charge exchange is typically nega- 
tive (repulsive, / < 0). However, the effective attractive 
of this form (/ > 0) is also possible [251 - 127] and in particu- 
lar it can originate from the coupling of electrons with in- 
tersite (intermolecular) vibrations via modulation of the 
hopping integral [25^ . or from the on-site hybridization 
term in generalized periodic Anderson model [261 [27] . 

The ferromagnetic XY order of pseudospins pi (for 
/ > 0) corresponds to the SS phase (s-pairing su- 
perconducting), whereas the antiferromagnetic XY or- 
der (for / < 0) - to the 77S phase (77-pairing super- 
conducting). However, in the absence of the exter- 
nal field in the XY-direction conjugated with the or- 
der parameter (A = A55 — ^ J^i {p7)- fo'' / > and 

A^s = F Sj ^^^P (**9 ■ ^i)iPi)- fo'' I <0, Q is a half of 
the smallest reciprocal lattice vector) there is a symmetry 
between the ferromagnetic and antiferromagnetic cases 
(except for an obvious redefinition of the order param- 
eter) for lattices consisting of two interpenetrating sub- 
lattices (such as for example SC or BCC lattices), thus 
we restrict ourselves to the / > case in the following. 
Let us stress that the symmetry is broken by the single 
electron hopping term i^ 7^ (cf . section VI I . 



The first analysis of the phase diagram of model (fll 
have been performed by Bari ^28j and Ho and Barry |.29] 
using the variational method in order to examine the in- 
stability of the Mott insulator to superconductivity for 
the special case of the half-filled band {n = 1). The ef- 
fects of diagonal disorder on the critical temperature for 
U — and n = I have been also determined [30', arriv- 
ing at a satisfactory qualitative interpretation of quite a 
number of different experiments in amorphous supercon- 
ductors. 

In the analysis we have adopted the variational ap- 
proach (VA) which treats the on-site interaction U ex- 
actly and the intersite interactions / within the mean- 
field approximation (MFA). We have also derived exact 
ground state results for a linear chain {d — 1). More- 
over, at T = several results obtained within the random 
phase approximation (RPA) and the spin-wave approxi- 
mation (SWA) (for d = 2 and d = 3 lattices) and within 
the systematic low density expansions (LDE) (for d = 3 
lattices) are presented. 

Within the VA the phase diagrams of model (fTl) have 
been investigated in [311 [52] as a function of the electron 
concentration n, but the stability conditions of states 
with phase separation have not been discussed. 

In this paper, we investigate the thermodynamic prop- 
erties of model Hamiltonian ([Tl for arbitrary chemical 
potential and arbitrary electron concentration. We fo- 
cus on the states with phase separation, which have not 
been analyzed till now. Arbitrary values of on-site inter- 
action U are considered and the range of pair hopping 
interaction / is restricted to the nearest neighbours. Our 
investigation of the general case finds that, depending on 
the values of the interaction parameters and the electron 
concentration, the system can exhibit the superconduct- 
ing ordered and nonordered homogeneous phases as well 
as the phase separation between them. Transitions be- 
tween various states and phases can be continuous and 
discontinuous, what implies existence of tricritical points 
on the phase diagrams. We present detailed results con- 
cerning the evolution of phase diagrams and thermody- 
namic properties as a function of interaction parameters, 
chemical potential p and electron concentration n. By 
taking into account the phase separated states and pre- 
senting several new results for lattices of various dimen- 
sions derived beyond VA (rigorous solution for a chain 
(d = 1), RPA for SQ and SC lattices and LDE in d = 3) 
our paper substantially extends and generalizes the re- 
sults and conclusions of previous papers concerning the 
model considered [31} [32] . 

The paper is organized as follows. In section |TT] we 
describe the method used in this work. There are also 
derived explicit formulas for the free energies of homoge- 
neous phases and states with phase separation as well as 
equations determining the superconducting order param- 
eter and the chemical potential in homogeneous phases 
obtained in VA. In section |III| we analyze the proper- 
ties of the system at zero temperature and present ex- 



act ground state results for a chain (section III A), dia 



grams obtained in the VA (section III B ) as well as re 



suits obtained in some approximations going beyond VA 
for d = 2,3 lattices (section IIIC and section HID). Sec- 
tion |IV| is devoted to the study of the finite temperature 
phase diagrams derived within the VA as a function of fi 
(section IV A ) and as a function of n (section IV B I . Some 
particular temperature dependencies of the thermody- 
namic parameters are discussed in section [V] Finally, 
section |VI| is devoted to conclusions and supplementary 
discussion. The appendices report some relevant details: 
explicit expressions of site-dependent self-consistent VA 
equations (IA]), particle-hole transformations (IbI, and se- 
lected analytical VA results (ICl). 



II. THE VARIATIONAL METHOD 

The free energy of the system as well as the self- 
consistent equations for the average number of electrons 
in the system and the average value of the order param- 
eter are derived within site-dependent VA in Appendix 
[K\ Restricting our analysis to the nearest neighbours and 
assuming no spatial variations of the order parameter the 
grand potential per site obtained in VA is given by: 

^A)-f = -A + 2/o|Ap-^ln(2Z), (7) 



where 



Z = cosh (^/3v/m^ + 4|/oA|2j + exp {I3U/2) , 
fl^^fi-U/2, Io = zI, A*-l^(p+), 



fin 



N 



and /? — l/(fcsT). The free energy per site: f ~ ut 
is derived as 

f{n) = fiin - 1) + {U/2)n + 2Io\A\^ - ^ ln(2Z). (8) 

The grand potential w is minimized in the system con- 
sidered for fixed fi, whereas the free energy / takes the 
minimum when n is fixed. 

The condition for electron concentration ([3]) and a min- 
imization of u! (or /) with respect to the superconducting 
order parameter |A| lead to the following self-consistent 
equations (for homogeneous phases): 



/isinh /37/i2+4|/oA| 



1 



sinh (Py/, 



4|/oA|- 



2+4|/nA|2 



Z^fl^+4\IoA\ 



0. 



(9) 



(10) 



In the homogeneous phases the double occupancy per site 



defined as D = -^ ^ 



D = 



fii-ffii^), has the following form 
^ lexp(/3[//2) 



(11) 



We also introduce the concentration of locally paired 
electrons Up = 2D and the ratio Up/n = 2D/n. No- 
tice that D is different from the condensate density 
(a fraction of pairs in the condensate) uq. These 
quantities differ even in the VA at T = 0: D ^ n/2 
and no = |(p^)P = (l/4)n(2 — n) in the superconducting 
phase. 

Equations ([9|-([To| are solved numerically for T > 
and we obtain | A | and n when /i is fixed or |A| and 
fi when n is fixed. The superconducting phase (SS) is 
characterized by non-zero value of |A|, whereas |A| = 
in the non-ordered (normal) phase (NO). 

Equation ( 10 1 is only the necessary condition for an ex- 
tremum of (|7| (or (|8|) thus the solutions of (l9|-(10l can 
correspond to a minimum or a maximum (or a point of 
infiection) of a; (or /). In addition, the number of min- 
ima can be larger than one, so it is very important to find 
the solution which corresponds to the global minimum of 

(0 (or (§). 

We say that the solution of (|9l-( 10 ) corresponds to 
a metastable phase if it corresponds to a minimum of 
the free energy / or grand potential w and the stability 
condition 



dfi 

dn 



>0 



(12) 



is fulfilled. Otherwise, we say that the phase is unstable. 
A stable (homogeneous) phase is a metastable phase with 
the lowest free energy (among all metastable phases and 
the phase separated states described below). 

Phase separation (PS) is a state in which two domains 
with different electron concentration: n+ and ri_ exist 
in the system (coexistence of two homogeneous phases) . 
The free energies of the PS states are calculated from the 
expression: 

fpsin+,n-) =mf+{n+) + (1 - TO)/_(n_), (13) 

where f± {n± ) are values of a free energy of two separat- 
ing phases at n± corresponding to the lowest homoge- 
neous solution for a given phase, m is a fraction of the 
system with a charge density n+ , 1 — rn is a fraction with 
density n_ (n_|_ > n_) and ttiti^ + (1 — m)n_ = n. The 
minimization of ( 13 1 with respect to n^ and ?i_ yields 



the equality between the chemical potentials in both do- 



mains: 



M+("-+) =M-("-) 



(14) 



(chemical equilibrium) and the following equation (so- 
called Maxwell's construction): 



li+{n+) 



/+(n+)-/_(n_) 



n+ — n- 



(15) 



which is equivalent with equality of grand potentials per 
site in domains: w-|_(/i_|_) ~ uj^{fi^). It implies that the 
transitions with a discontinuous change of n in the sys- 
tem considered for fixed fi can lead to occurrence of the 



regions of phase separation in the concentration range 
n_ < n < n_|_ on the diagrams obtained as a function 
of n. In these regions the homogeneous phases can be 
metastable as well as unstable, depending on the n- 
dependence of /i (cf. (12l). In the PS states the chem- 
ical potential // = /i_|_(n+) = /i_(n_) is independent of 
the electron concentration, i.e. dfi/dn = 0. For more de- 
tailed discussion of Maxwell's construction we refer the 
reader to existing literature, e. g. [Ml[5?| . 

In the model considered only one type of PS states can 
occur, which is a coexistence of SS and NO phases. 

In the paper we have used the following conven- 
tion. A second (first) order transition is a transition 
between homogeneous phases with a (dis-)continuous 
change of the order parameter at the transition temper- 
ature. A transition between a homogeneous phase and 
the PS state is symbolically named as a "third order" 
transition. At this transition a size of one domain in 
the PS state decreases continuously to zero at the transi- 
tion temperature. Second order transitions are denoted 
by solid lines on phase diagrams, whereas dotted and 
dashed curves denote first order and "third order" tran- 
sitions, respectively. 

The phase diagrams obtained are symmetric with re- 
spect to half-filling because of the particle-hole symme- 
try of Hamiltonian (IT|) , so the diagrams will be presented 
only in the range // < and < n < 1. The transforma- 
tion formulas of various thermodynamic parameters and 
quantities such as electron concentration, superconduct- 
ing order parameter, chemical potential and free energy 
under the electron-hole transformation are collected in 
Appendix |B] 

The numerical computations have been done in 
C/C++ with source codes from J3S|- The minima have 
been found by using Brent's method in one dimension 
and downhill simplex method in multidimensions. We 
have also implemented Newton-Raphson method using 
derivative to solve non-linear equations. 



III. THE GROUND STATE 

Let us consider the ground state (GS) of model (fTl . 

In the absence of tij term, the parity of the number of 
particles at each site is a good quantum number. Thus, 
one can decompose the eigenspace of H into sectors spec- 
ified by the parity of the occupation number at each site 
(n.i). Within any sector, the system is a collection of even 
segments, within which pairs are confined (the subspace 
with excluded singly occupied sites) interleaved by odd 
segments (the subspace with excluded doubly occupied 
sites) . 

Since 1^0, the interfaces between the odd and even 
regions will cost energy, so in the thermodynamic limit 
the ground state will be either a single infinitely large 
even segment or a single odd segment for fixed /i. Only 
in the case when the grand potentials per site (it is the 
chemical potential) of the segments are equal there is 



possibility of coexistence of both types of segments if 
N — >■ +CX) - it appears as a first order line on the GS 
phase diagram. 

The situation is different when n is fixed. In such case 
the PS states can occur in definite ranges of n. The rea- 
son for this can be explained as the following. Let the 
energy of a single boundary (interface) between odd and 
even segments be Ej. If the number of intersite bound- 
aries will increase as N'' (where 7 < 1) the contribution 
from them to the total free energy N^Ej/N — Ei/N^~^ 
disappears in the thermodynamic limit N — > +00. 

Within the subspace excluding the double occupancy 
of sites (by electrons for n < 1 or holes for n > 1) the 
product pf pj ~ for any i ^ j, i. e. the / term does 
not contribute. Effective Hamiltonian (fTl) for this case 
takes the form Hj = U ^^ fii^fii^ — p, ^^ fii and the cal- 
culations of the GS grand potential to mo — {Hi) /N and 
the GS energy E^o = {Hi + pJ2i "-i)/^ ^re quite obvi- 
ous. As a result in the limit /? — )■ +00 from (IS]) and (l9| 
one obtains for fixed n in the NO phase: 



ENo{n) =0, for n < 1 

ENo{n) = U{n — 1), for n> 1 
Eno ~ 0, for n = 1 



(m = -C//2), 

(m = C//2), (16) 

(m = o). 



This GS is infinitely degenerate for N ^ 00. For n < 1 
(n > 1), nN ((2 — n)N) sites are occupied by electrons 
(holes) and the spins of all particles are independent 
of each other. Inclusion of the intersite spin-spin and 
density-density interactions can remove this degeneracy 
and generate various magnetic and charge orderings. The 
magnetic moments are not suppressed in this state. For 
ri 7^ 1 there is no order pattern in the distribution of 
the particles, whereas for n = 1 Mott phase with one- 
electron-per-site configuration (with random direction of 
electron spin per site) occurs. 

From (l7| and ^ for fixed p, and /3 -^• +00 one obtains: 



U)' 



NO 

"^NO 



•"NO 



{p) =0 (n = 0, NO - empty), 

Ip) = ~p-U/2 (n = 1, NO - Mott), (17) 

{p) = -2p (n = 2, NO -full). 



The average occupancy of sites can only be integer for 
fixed p, what can be intuitively understand if one re- 
alizes that {fii) = 0,1,2 only and every site is equiva- 
lent (there is no interactions between different sites in 
Hi). The state w^q ^^ ~ ^™Pty) is stable with respect 
'^o (NO - Mott) if [/ < -2/2 (and /2 < 0), whereas 



to w 



^NO 



< uj%Q liU <2p (and p > 0). 



Within the subspace excluding the single occupancy of 
sites the eigenvalues of charge operators are pf = ±1/2, 
(Pi)^ — 1/4 for each site. Thus, from (W]), the effective 
Hamiltonian for this case takes the form: 



Hii^~lY.iPtPj+PtPi) 



mE(2p? 



1) 



(18) 



(iJ) 



with auxiliary condition ([6| . Hn has the form of the XY 
quantum spin model (S = 1/2) with an effective external 



field Jl = fji — C//2 in the z-direction, such that the aver- 
age magnetization has a fixed value equal to [n — l)/2 
(cf. (|6])). The GS energy of Hamiltonian (18 1 is derived 
as 



1 



1. 



Ess{.n) ^ -{Hi I + ^iY^fi,) = Exvin) + -Un, (19) 



where ExY{n) = -jf{J2{^J) {ptPj + P'^ P^ ))■ 
grand potential in this case has the form 

^ss{p) = ExY{n{fl)) - p,n{p,), 



The GS 



(20) 



where it is underlined that /l is an independent variable. 

As it follows from previous analyzes [TT, T^i the GS of 

model ( 18 ) exhibits the long-range XY order for d> 2 



dimensional lattice in the whole concentration range 
< n < 2. For d — I, only the short range XY correla- 
tions are developed [SB] , nevertheless even arbitrary small 
interchain interactions can induce the real long-range or- 
der also in this case. 

The SS phase will be realized if E^o > Ess for fixed n 
(if wtvo > ^ss for fixed p.). For explicit calculation of the 
phase boundary between the NO (Eqs. ( 16 ) and ( 17 1) and 



the SS (Eqs. (19 1 and (pO|)) phases and determination of 



the GS properties of the SS phase for various dimensions 
and the electron concentration or the chemical potential 
one has to find the solution of model ( 18 1 at T = Q. 

Let us stress that if n is fixed one has to consider also 
the PS state in which one domain is in the NO phase 
and the other is in the SS phase and which energy is 
calculated from (13 1. However, the GS energy of the 



PS (SS-NO) state can be obtained in the most cases 
only numerically by minimization of the free energy given 
by ( [T3| . In the NO domain njvo = 1 (Mott phase with 
one-electron-per-site configuration), whereas concentra- 
tion in the SS domain is dependent on U/Iq (and ap- 
proach used). 



Exact results in d - 



and the exact expressions for chemical potential ^ and 
the superfiuid density j^ (the short-range order parame- 
ter) are: 



P'SS 



dEE 



J 



dn 

v = mp. 



^TT 17-1 '^'^ 

- i^U = -\Io\cos—, 



i+ii 



1 7rn 

TT Z 



(23) 

(24) 



where i + 1 denotes the nearest neighbour of the i-site. 
For /o > (SS case) ferromagnetic correlations between 
XY components of neighbouring pseudospins are present 
{Pi Pi+i) > 0' whereas for /q < (r/S case) antiferromag- 



netic correlations occur {pfp~^i) < 0. From (23 1 one ob- 
tains that the SS can exist in the range —1 <p^\Io\ < 1. 
Notice that stability condition ( [l2| is always fulfilled in 
the SS phase. 

The GS phase diagrams obtained are shown in figure^ 
as a function of fl (panel (a)), and as a function of n 
without (panel (b)) and with (panel (c)) consideration of 
the PS state. At half-filling there is no phase separation 
and the NO - Mott phase is stable for U/Iq > 4/7r. The 
boundaries between different states corresponding to the 
exact results in d = 1 are denoted by dotted lines (if they 
are different from those obtained in VA) . Notice that for 
homogeneous phases dfi/dn > in the ranges of the PS 
state occurrence. 



B. Variational approximation (d -^ +oo) 

For 1 < d < oo there are no exact solutions for model 
([l8| av ailable (except of the LDE for d — i - see |13l 
I14 [ I31 | ) and the approximate approaches are necessary 
to find Ess- The simplest of these approaches is the VA 
which at T = gives (the /3 -^ +oo Hmit of ([8|-([l0| for 
A^O) 



Ess^\un-]^\h\n{2-n), 
|A|2=j, = ^j(2-n), 



(25) 
(26) 



The exact solutions of model ( 18 1 for arbitrary n can 
be derived for d = 1 and d = -|-(xi. For d = 1, by making 
use the exact results for the ground state energy of d = 1 
XY model in a transverse field fSSHSSl one obtains: 



ExY - -2|/o 



1 



— sm — 

TT 2 



(21) 



where /q — 21 (because z = 2). Notice, that due to a well 
known isomorphism between the planar ferromagnet and 
antiferromagnet for loose packed (alternating) lattices 
with nearest-neighbour |40l HI] , the results for Ejfi^ are 



the same for both signs of I. Thus, from (19 1 and (21) 
we have 



Eh'^ = \un 



^1 T , 1 . Trn 
2|/o|-sm — 

TT Z 



(22) 



and 



fJ-ss 



-|/o|(l-n). 



(27) 



Equation ( |27[ ) implies that the SS phase can exist in the 
range —1 < p-/\Io\ < 1. Notice that Up = 2D — n, which 
means that all electrons in the system are locally paired 
in the SS phase at T = 0, whereas the finite value of the 
depletion n' — n — 2no {uq = |Ap) at T = is due to the 
hard-core effect |T^ . 

The GS energy of the PS state can be obtained by 
minimization of ( 13 1. In the NO domain nxo = 1 (Mott 



phase), whereas concentration in the SS domain is depen- 
dent on U/Io and nss = l± V^/Io - 1 (1 < U/Iq < 2). 
It corresponds to the first order SS-NO boundary on 
the U/Iq vs. fi/Io GS diagram determined by equation 
(/2//n)2 + 1 = U/Io. 
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FIG. 1. (Color online) Ground state phase diagrams: (a) as a function of p,/Io, (b) and (c) as a function of n without and with 
consideration of the PS states, respectively. Dotted, dashed, dashed-dotted and solid lines denote the boundaries for d = 1 
(exact results), d — 2 (RPA, SQ lattice), d = 3 (RPA, SG lattice), and d — > +00 (VA result), respectively. At half-filling (n = 1) 
the NO (Mott state) is stable above the end of the PS-SS line. The NO-SS transition for p./Io ~ —i (panel (a)) is second 
order, all other transitions between homogeneous phases (panels (a), (b)) are first order. 



The GS phase diagrams obtained in VA are shown 
in figure ^ The VA phase boundaries are denoted by 
soHd Hnes. Similarly as in section |III A| one finds that 
d(x/dn > for homogeneous phases in the ranges of the 
PS state occurrence. 

For hypercubic lattice of the dimension d = 00 the VA 
becomes the exact theory for model (IT|) (as well as for 
model ([18|). 



C. Random phase approximation {d — 2,3) 

Going beyond VA for 1 < d < +00 a reliable approach 
is the self-consistent RPA [T3t [M] , This approach has 
been proven to be a very good approximation scheme 
in problems of quantum magnetism and it fully takes 
into account quantum fluctuations, which can be of cru- 
cial importance for the considered system for d < 3. 
The GS energy of model ([l8|) is given by (19 1, where 
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EkhiPkPk: 

space-Fourier transform of the charge operators and the 
pair hopping integral, respectively, and the chemical po- 
tential is determined by (n — l)/2 = (1/-/V) J^iiPf) 

The energy E^y^ is obtained by means of Green's 
functions using the spectral theorem, and by applying 
the approximation for the longitudinal correlation func- 
tions [II]. At T = this yields 
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the self-consistent equation 

^"' = ^E^ = 1 + 2^0' (29) 
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(30) 



The result corresponding to the SWA (being the lowest 
order expansion of above equation) is i? = 1 — 2ipQ with 
2V'o determined by ( 30 1 . 



For the SS phase the spectrum Ej^ is gapless and linear 

in |fc| for small k. One gets E^^^ = \k\s with the sound 
velocity s given by s = 2RIay/z sin 9 {h = 1). 

The chemical potential in the SS phase in the RPA is 
given by p.§g = Io{n — 1), which is the same as in VA. 
The quantum corrections to this VA result for p,ss can 
be evaluated as in [T3] and one obtains: 
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where e^ = y^(l - 7^)' + 7^(1 - 7*:) sin^ 0, 7^ = /g//o. 

One should notice that RPA and SWA give very similar 
result for the model ^ in the GS [HIIS]. The result- 
ing phase boundaries for d = 2 square (SQ) lattice and 
d = 3 simple cubic (SC) lattice are shown in figure [l] and 
denoted by dashed and dashed-dotted lines, respectively. 



D. Low-density expansion for d = 3 lattices 



^k = H 



%- 



Bj: = RI^ sin^ 



where 

Ej: ~ \i ■^- ^ B?. is the collective excitation spec- 
trum, eS — 2(/o — Ij:), is given by cos^ 6 = {n— lY /R^ 
and the length of pseudospin R is given as a solution of 



In |13l [H] several rigorous results concerning the GS 
properties of the model of hard-core charged bosons have 
been derived for various lattices in d = 3 using a system- 
atic LDE based on the knowledge of the exact two-body 
scattering amplitude. Transforming the results of sec- 
tion VI B from [14 into our problem one gets for n — >■ 




Ti/I 








FIG. 2. (Color online) The SS-NO transition temperature 
as a function of p./ lo for increasing values of U / lo (numbers 
above the curves). Solid and dotted lines indicate second 
order and first order transitions, respectively. The dashed- 
dotted curve (labeled as TCP) is a projection of the tricritical 
point line on the ksT/Io vs. fl/Io plane. 



FIG. 3. (Color online) Finite temperature phase diagram 
plotted as a function of fi/Io and U/Iq. Solid and dotted 
lines indicate second order and first order transitions, respec- 
tively. The tricritical point line is denoted by dashed-dotted 



the following exact expressions. The GS energy of SS 
phase Egg^ is given by 
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where a = 1/C is the exact scattering length for the con- 
sidered case, C = (l/N) J2ki'^ - ^fe/^o)"^ is the Watson 
integral for a given lattice, 7 — Vo/a'^, Vq is the vol- 
ume of the unit cell, a is the lattice constant. In par- 
ticular, for SC: C = 1.51638, z = 6, 7 = 1; for BCC: 
C = 1.3932, z = 8, 7 = 1/2; and for FCC: C = 1.3446, 

z =: 12, 7 = 1/4. 

The chemical potential has a form 
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E. Summary of the GS results 

The general structure of the GS phase diagrams de- 
rived is similar for all lattice dimensionalities (d > 1) and 
approaches used. There are two sequences of phase tran- 
sitions with increasing U/Iq for fixed n: 

(i) SS^PS^NO (forn^l), 

(ii) SS -^ NO (Mott) (for n = l). 

The effects of quantum fluctuations increase with de- 
creasing lattice dimensionality and for a linear chain their 
influence on phase diagram is the most prominent. Thus 
the region of the homogenous SS phase stability decreases 
with increasing d. In d — > +cx} the quantum fluctuations 
are suppressed. 

One should also notice that the transition PS-NO at 
U/Iq = 2 is independent of n and dimensionality of the 
lattice. The PS region expands with increasing dimen- 
sionality. 



and the sound velocity s is derived by 



fi^s^ ^UoIa^'n{2-n)a 
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(34) 



The phase diagrams obtained within LDE are in very 
good agreement with RPA results (for SC lattice) in the 
range n < 0.2. The homogeneous SS phase boundary de- 
rived within LDE is moved little downwards in compari- 
son with the RPA ones. An analysis of LDE results ob- 
tained for various cubic lattices in d = 3 shows that the 
region of the SS phase stability shrinks with increasing 
number of the nearest neighbours in the lattice. 



IV. FINITE TEMPERATURE PHASE 
DIAGRAMS DERIVED WITH VA 

A. Analysis for fixed /j, 

All phase transition boundaries, necessary to construct 
the complete phase diagram for fixed /i, have been ob- 
tained numerically by solving the set of equations ^- 
( 10 1 and comparing the grand potential ^ for the so- 
lutions found. The resulting ksT/Io vs. pjla phase di- 
agrams obtained for various fixed values of the on-site 
interaction U/Iq are shown in figure p] The overall be- 
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FIG. 4. (Color online) Phase diagrams fesT/Zo vs. n for different values of U/Iq (as labeled). Solid and dashed lines indicate 
second order and "third order" boundaries, respectively. T denotes the tricritical point. 



haviour of the system as a function of ^ is demonstrated 
in figure [3] 

For U/Iq < (2/3) In 2 only the second order SS-NO 
transitions occur with increasing temperature. The tran- 
sition temperature is maximal for U — > — oo, p, = and it 
decreases monotonically with increasing U/Iq and \p,\/Io. 

The most interesting is the range 
(2/3) In 2 < U/Iq < 2. In this range the /i-dependence 
of the SS-NO transition temperature becomes non- 
monotonic and with increasing U/Iq the maximum of 
the transition temperature moves towards higher values 
of |/i|//o- In this range of the on-site interaction, the 
order of the SS-NO transition changes from second 
order to first order with increasing of fi/Io- A tricritical 
point (TCP) is connected with the change of transition 
order. 

For 2 < U/Iq < +cx) only the NO phase is stable at 
any T > 0. 

As it will be showed in section IIVBI the first order 
phase boundaries are associated with the existence of the 
PS states in definite ranges of n if the system is consid- 
ered for fixed n. 



B. Analysis as a function of n 

The transition temperatures have been determined by 
self-consistent solving of (|9|- ([To| and ([l4|-( 15 1 and com- 
paring free energies 



and ( |13p . 

The concentration dependencies of (i) the second or- 
der SS-NO transition temperatures and (ii) the "third 
order" SS-PS transition temperatures can be also easily 
obtained from the diagrams as a function of p,, because 



of a simple relation pss — Io{n ~ 1) between the electron 
concentration and the chemical potential in SS phase. 
The transitions (i) and (ii) for fixed n correspond to the 
second order and first order SS-NO transitions for fixed 
/i, respectively. However, there is no simple correspon- 
dence for the PS-NO transitions (cf. (|CT|). 

Examples of the fc^T/io vs. n phase diagrams evalu- 
ated for various values of the on-site interaction U/Iq are 
shown in figure |4] 

For U/Iq < (2/3) In 2 the PS states do not occur and 
the obtained phase diagrams are the same as those de- 
rived in |31] . The transition between homogeneous SS 
and NO phases taking place with increasing temperature 
is second order for arbitrary n. The SS-NO transition 
temperature is maximal for [/ — > — cx), ti = 1 and it de- 
creases monotonically with increasing U/Iq and decreas- 
ing n (for n < 1) - cf. figure l4k. 

In the range (2/3) ln2 < C//Jo < 2 with increasing 
U/Iq the maximum of the SS-NO transition tempera- 
ture moves towards lower concentrations (for n < 1). In 
a definite range of U/Iq and n the PS state SS-NO is 
stable. For 1 < U/Iq < 2 the PS state extends from the 
ground state (cf. figures |4li, e, f and figure fl|:), whereas 
for (2/3) In 2 < U/Iq < 1 it is stable only at finite tem- 
peratures (cf. figures Ep, c). The critical point for the 
phase separation (denoted as T, which is a tricritical 
point, TCP) lies on the end of the SS-NO second order 
line. 

Between U/Iq = 2 and U/Iq — > +oo only the NO phase 
is stable at any T > and there are no transitions. 

For given value of n the UbT/Iq vs. U/Iq phase dia- 
grams of the system are shown in figure [Sj All of the 
transition temperatures decrease with increasing U/Iq. 
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FIG. 5. (Color online) Phase diagrams ksT/lQ vs. f///o for n = 0.25 (a), n = 0.75 (b), and n= 1 (c). Solid and dashed lines 
indicate second order and "third order" boundaries, respectively. T denotes the tricritical point. 
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FIG. 6. (Color online) Temperature dependencies of (a) the superconducting order parameter |A| (b) the ratio np/n and 
(c) the specific heat c/fcs plotted for: Ji/Iq = -0.1 and U/Io = -0.5 (A), U/h = -0.25 (B), U/h = 0.75 (C) as well as for: 
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In the range < n < 1 the diagrams are rather similar 
and they consist of three regions (exemplary diagrams are 
shown in figures pa, b). The regions of SS and NO phases 
stability are separated by second order boundary or, for 
larger on-site repulsion, by the region of PS state occur- 
rence. For half-filling {n — 1, ft = 0) only homogeneous 
phases are present on the diagram (figure Isb) . The T- 
point, connected with the change of the SS-NO transition 
order, is located at ksT/Ia = 1/3 and U/Ia = (2/3) In 2 
(for n — 1). 

The possible sequences of transitions with increasing 
temperatures and the transition orders of them are listed 
below: 

(i) SS^-NO: second order for n 7^ 1 and second order 
or first order for n — 1 (and U/Iq < 2), 

(ii) PS^-NO: "third order", it can take place only for 
n^ 1 (and 1 < U/Iq < 2), 

(iii) SS— >-PS— ^NO: both "third order", it can take place 
only for n ^ 1 (and (2/3) ln2 < U/Iq < 2). 

Notice that increasing temperature usually suppresses 
heterogeneity in a system what favours occurrence of ho- 
mogeneous phases rather than PS states. In the model 
considered the PS state can develop at higher tempera- 
tures than the homogeneous SS phase. This rather un- 
usual behaviour can be connected with the interplay be- 



tween the pair hopping interaction and the on-site repul- 



V. THERMODYNAMIC PROPERTIES (VA) 

In this section we present several representative depen- 
dencies of the thermodynamic characteristics for fixed 
model parameters (figures [6}|8| . 

In particular for fixed jl (figure pi, one can single 
out two limiting types of thermodynamic behaviour near 
transition temperature Tss'- (i) the local pair regime 
and (ii) the pair breaking regime. In between, there is 
a crossover between the two regimes (see e. g. the plots 
rip/n vs. kBT/Io in figure pp). 

For large on-site attraction the concentration of locally 
paired electrons Up = 2D exhibits no sharp feature as the 
temperature is lowered through Tss- The number of non- 
paired electrons at Tss is negligible and the transition 
is to the state of dynamically disordered pairs (only for 
\U\/Io 3> 1 and [/ < 0, the local pair regime). 

In the second limit, for on-site repulsion U < 2/o, rip 
has a sharp break at Tss and a substantial fraction of sin- 
gle particles can exists above Tss- We call this the pair 
breaking regime. As temperature is lowered, the conden- 
sate growths both from a condensation of pre-existing 
pairs and from binding and condensation of single parti- 
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FIG. 7. (Color online) Temperature dependencies of (a) the specific heat c/ks, (b) the ratio rip/n, and (c) the fractions 
m-i and m2 of the system with a charge density n_|. and n_, respectively. On the insets there are temperature depen- 
dencies of the entropy s/fcs (a) and the chemical potential fi/Io (b). The characteristics are obtained for: n = 0.8 and 
U/h = 1.25 (A), n = 0.7 and U/Io = 1 (B), and n = 0.2 and U/h = 0.75 (C). 
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cles. For small binding energies, ifn<^ 1 (/x « —1), there 
will be essentially no pre-formed pairs at Tss- 

Notice that non-zero value of Up does not imply that 
local pairs are in coherent state and even significant val- 
ues of Up/n are possible in the NO phase. In the limit 
T -> -l-oo Up increases to Up/n — >■ 0.5 (each of four states 
at a given site can be occupied with equal probabil- 
ity). The condensate density (which can be approxi- 
mated as no « |Ap at least for n <^ 1, uq y^ D) vanishes 
for T > Tss, but the doubly occupied states are still ther- 
mally excited above Tss {D 7^ 0). 

In figure [6^ the temperature dependencies of the su- 
perconducting order parameter A are presented, where 
one can see clearly the discontinuous change of the order 
parameter (lines C and G). The other lines correspond 
to second order transitions. 

Finally, let us briefly summarize the behaviour of the 
speciflc heat at constant volume c = — T ( 4^^ J (flg- 

urelot). The NO phase is characterized by the relatively 
broad maximum in c connected with continuous changes 
in a short-range electronic ordering (in higher tempera- 
tures, not shown in figure^. The narrow peak in c{T) 
is associated with the first order transition, while the A- 
point behaviour is typical for the second order transition. 
Let us concentrate now on the thermodynamic 



properties of the system in the PS state and their 
changes at the "third order" transitions for fixed elec- 
tron concentration n (figures [t] and jsl. As an ex- 
ample, we consider the temperature dependencies of: 

(i) the specific heat c — —T ( ^^4 ) , (ii) the entropy 

s — —-Jp, (iii) the concentration of paired electrons np/n, 
(iv) the chemical potential /i//o, and (v) the fractions 
TOi — {n — n^)/{n+ — n_) and 7712 = 1 — mi (sizes of do- 
mains) with a charge density n^ and n_ (figure^ as well 
as the concentration dependencies of: (a) /i//o, (b) f/Ia 
and (c) Up/n (figure [8|. 

Thermodynamic characteristics for three different pos- 
sible sequences of transitions with increasing tempera- 
ture are shown in figure It] The label (A) corresponds to 
the case of PS— >N0 transition. At such a transition the 
specific heat c exhibits a finite jump at transition tem- 
perature, what is associated with the continuity of the 
entropy s (cf. figure l7k). The lines labelled by (B) cor- 
respond to SS^-PS-^NO case. At the SS-PS transition 
there is also a finite jump in c, however the c in the PS 
state is larger than that in the SS phase, contrary to the 
PS-NO transition, where the jump has an opposite sign. 
Finally, the lines (C) correspond to the case of second 
order SS^-NO transition. The concentration ratio of lo- 
cally paired electrons Up is presented in figure l7b. On 
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the inset of figure [7]d the chemical potential is shown as 
a function of temperature. One can see that /2 in the SS 
phase is independent of T (cf. (C2|). In figure l7t there 
are shown dependencies of the fractions of the system 
TTii and 7712 with charge density n+ and n_, respectively. 
For instance, in the case (A) the SS domain with lower 
concentration (n_) vanishes at "third order" transition. 

The chemical potential fl/Io and the free energy 
f/Io as a function of concentration for U/Iq = 1.25 and 
ksT /Iq = 0.05 are presented in figure [sk and figure Isb, 
respectively. For n < n^ the SS phase is stable and its 
free energy is the lowest one. In the range n_ < n < n^- 
the PS state has the lowest energy (solid line). The so- 
lutions of (l9|-(10l corresponding to the lowest energy 



ones (the homogeneous phases with the lowest energy) 
are thermodynamically stable i. e. dp,/dn > in their ex- 
istence ranges (denoted by dashed lines in figures pa, b) . 
At point labelled as X the transition between metastable 
phases occurs, eg. for n < nx the free energy fss of 
the SS phase is lower than the energy Jmo of the NO 
phase, whereas fss > Jno for n > nx- The dotted lines 
in figures [8^, b correspond to the homogeneous phases 
with the highest energy. One should notice that in the 
described case for U/Iq — 1.25 and ksT/Io — 0.05 in the 
whole range of the PS state occurrence two homogeneous 
phases are metastable, whereas in the ranges of the ho- 
mogeneous phases stability there is one metastable phase 
at the most. The linear dependence fl/Io = n — 1 is char- 
acteristic for the SS phase, whereas in the PS state j2/Io 
is independent of n. This is consistent with square depen- 
dence of fss and linear dependence of fps as a function 
of n. 

The anomalous decrease of the chemical potential with 
doping in the lowest energy metastable homogeneous 
phases near X-point (i. e. an abrupt drop at X-point) 
is a signal that it is thermodynamically more convenient 
to phase separate into two subsystems with the different 
electron densities n_ and n+ (figure ^). The density 
range for which the PS occurs can be readily found by 
the Maxwell's construction described in section HIl 

The n-dependencies of the chemical potential /l//o, 
the free energy ///q and the concentration of paired 
electrons Up/n — 2D/n are shown on the insets of fig- 
ure [8^ and figure [8]3 and in figure [8J3, respectively, for 
U/Iq = 1.25 at several fixed temperatures. They are 
obtained for: keT/Io = 0.025 (A), fesT/Jo = 0.05 (B), 
ksT/Io = 0.075 (C), and keT/lQ ^ 0.1 (D). 

At "third order" transitions (SS-PS and PS-NO) the 
thermodynamic parameters such as fl, rip/n, s, / are 
continuous, whereas the c exhibits a finite jump as at 
the second order transitions. However, one should notice 
that the order parameter for "third order" transitions is 
the concentration difference n-|_ — 7i_ (not A, which is the 
order parameter in one domain) and its change is discon- 
tinuous at transition temperature. Such transitions are 
present if the system is consider for fixed n and they are 
associated with first order transitions at fixed /2. 



VI. FINAL REMARKS 

We have studied a simple model of a superconductor 
with very short coherence length (i. e. with the pair size 
being of the order of the radius of an effective lattice 
site) and considered the situation where the single parti- 
cle mobility is much smaller than the pair mobility and 
can be neglected. We have evaluated the phase diagrams 
and thermodynamic characteristics of the model at fixed 
/i and at fixed n and determined the ranges of the PS 
occurrence in the system. 

One should stress that within the VA the on-site U 
term is treated exactly. Thus, the major conclusions of 
our paper concerning the evolution of the properties of 
the system with U are reliable for arbitrary U. The mean- 
field used for the intersite term is best justified if the lij 
interactions are long-ranged or if the number of nearest 
neighbours is relatively large. The derived VA results 
are exact in the limit of infinite dimensions d — >■ -l-oo, 
where the MFA treatment of the intersite interaction / 
term becomes the rigorous one. Moreover, the VA yields 
exact results (in the thermodynamic limit) for /y of infi- 
nite range {lij = {l/N)I for any {i,j)) regardless of the 
dimensionality of the system [29J . 

The properties of the system are strongly dependent on 
the ratio of U/I, which is related to the relative values of 
the pair binding energy Eb and the effective pair mobil- 
ity tp. For the case of a single pair in a system, i. e. for 
n — > 0, we have Ei, = —U + 2/o, tp ~ 21 (the pair effec- 
tive mass is m* — \/{2tp) = 1/(4/)) [31] and increasing 
n can modify E}, . There are two well defined limits of the 
model: (i) the local pair regime and (ii) the pair breaking 
limit. In between there is a crossover between the two 
regimes for fixed ji. 

The local pair limit is realized for large pair binding 
energy (i. e. for C/ < 0). In such a case, the two inter- 
actions: the on-site interaction U and the intersite pair 
hopping / cooperate and transition temperature is deter- 
mined by pair mobility (i. e. by /). The SS-NO transition 
with increasing temperature is second order and is to the 
NO phase being a state of dynamically disordered local 
pairs. 

The opposite regime i. e. the pair breaking limit is re- 
alized for substantial values of f/ > 0. Repulsive U desta- 
bilizes on-site pairing and competes with /. As we have 
shown above, the increasing U: (i) for fixed /i it changes 
first the nature of phase transitions from a continuous to 
a discontinuous type and then it suppresses superconduc- 
tivity for |/i|//o = 1*^ ~ 1| ~^ 1; (ii) it stabilizes the phase 
separation state (SS-NO) in definite ranges of n and tem- 
perature. This behaviour is associated with the presence 
of the tricritical point on the phase diagrams. With in- 
creasing U/Iq the maximum of the transition tempera- 
ture moves towards lower concentrations (for n < 1). In 
definite range of U/Iq and n the PS state SS-NO is stable. 
For 1 < U/Iq < 2 the PS state extends from the ground 
state, whereas for (2/3) In 2 < U/Iq < 1 it is stable only 
at finite temperatures. 
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Finally, for large U {U/Iq > 2), the system remains 
in a normal state at all temperatures for any chemical 
potential or concentration. 

Comparing the VA results at T = with the other ones 
(exact ones for d = 1 chain, RPA ones for d = 2 square 
lattice and d — 3 cubic lattices) , one can see that the 
general structure of the diagrams is similar. It implies 
that the VA can give qualitatively reasonable results also 
for lattices of finite dimensionality. However, one should 
stress that in finite dimensions due to quantum fiuctua- 
tions connected with lij interactions the regions of the 
homogeneous SS phase occurrence are extended in com- 
parison with the VA results. Moreover, in a linear chain 
only the short-range order occurs at T = in contrary to 
the VA and RPA (SWA) in d > 2, where the long-range 
order is present. 

One should emphasize that for d < -l-oo the VA overes- 
timates transition temperatures and gives incorrect con- 
centration depe ndence of them {kBTss ~ — 1/ ln(n/2) for 



n ^ 1, cf. (C3l). For d = 3 ideal Bose gas with k spec- 



trum (the [/ ^ and n <^ 1 limit of the model consid- 
ered) the n-dependence of the SS-NO temperature is 
Tss ^ n^'^, what can be obtained, for example in the 
RPA, even for finite values of U/Iq [H i^. The RPA 
predicts also the existence of tricritical point in finite 
temperatures |32] in d = 3, what confirms that the VA 
can give qualitatively reasonable results also for real lat- 
tices in finite temperatures. The RPA results for d = 2 at 
T > [131 112] are in agreement with the Mermin- Wagner 
theorem |33] , which proofs the absence of the long-range 
order for d < 2 at T > 0. In d = 2 one has the Kosterlitz- 
Thouless-type transition with increasing temperature. 

The model considered can describe the systems in 
which the single electron mobility (single electron hop- 
ping) is much smaller than the pair mobility (kinetic en- 
ergy of pairs) and it can be neglected in first approxima- 
tion. Let us comment on the finite electron bandwidth 
effects. The presence of the hopping term "^^ j „ ^ij^ta'^jcr 
breaks a symmetry between / > (favouring SS) and 
/ < (favouring rjS) cases. The phase diagrams of the 
Penson-Kolb-Hubbard model are more complicated than 
those obtained for Hamiltonian (fTl) |16ff23| . even in the 
ground state. We can suppose that small but finite single 
electron hopping tij will not qualitatively alter the phase 
diagrams, at least for the case ksTss > X^j ^Ij/U- The 
main effect oitij (for |t/| ^ i^j, C/ < 0) is a renormaliza- 



W/|I| 



CO 





/ 1.5 2.0 

/ NO 



(Mott) 



PS [NO/NO] 



tion of the pair hopping term lij -^ lij 
introduction of an effective intersite density-density re- 
pulsion ^ tfJ\U\. For U < and / < the charge den- 
sity wave state can also occur |171 FIB] . For U > and 
both signs of / the tij term generates magnetic correla- 
tions competing with superconducting ones and its effects 
essentially modify the phase diagrams and the proper- 
ties of normal state (particularly for larger values of tij). 
In such a case it is necessary to consider also various 
magnetic orderings. Moreover, several phase separation 
states involving superconducting, charge and (or) mag- 
netic orderings could also be stable for n ^ 1. For larger 



FIG. 9. Ground state phase diagram in presence of density- 
density interaction W between nearest neighbours for half- 
filling obtained within VA (solid lines). The dashed lines de- 
note the rigorous boundaries for ID-chain. All transitions 
between homogeneous phases are discontinuous. 7o = zl. 



values of tij the studies of GS phase diagrams of the half- 
filled PKH model for d-dimensional lattices have been 
determined by means of the (broken symmetry) Hartree- 
Fock approximation (HFA) and by the slave-boson mean- 
field method (SBMFA) in [T7]. For d = 1 the diagrams 
supported also by the results of the continuum-limit field 
theory and exact diagonalisation studies [18, are shown to 
consist of at least nine different phases including super- 
conducting states, bond located antiferromagnetic and 
charge-density-wave states as well as mixed phases with 
coexisting site and bond orderings. The stability range 
of the bond-type orderings shrinks with increasing lat- 
tice dimensionality and for d — > -|-oo the phase diagrams 
involve exclusively site-located orderings. One should no- 
tice that the phase boundary at T = and n = 1 deter- 
mined by the HFA and SBMFA {tij ^ 0) are in agree- 
ment with our VA result U/Iq = 1, even for relatively 
large tij (especially for / > 0) [17] . Obviously, for tij ^ 
the transition is to the antiferromagnetic Mott phase. 

The intersite density-density interactions Wij can in- 
troduce mixed phases and (or) phase separations between 
charge orderings and superconductivity [44] . The ground 
state phase diagram in presence of interaction W between 

■hihi to 



2i^ /IC/I and an nearest neighbours (added a term {W/2)Y^ 



Hamiltonian (Jlj)) calculated for half-filling is presented 
in figure [9] (cf. with [15], where the case tij ^ has 
been analyzed). For sufficiently strong intersite repulsion 
M^ > the homogeneous charge ordered (CO) phase oc- 
curs (with uq = (1/-/V) X^ii^-i) 6^P (*Q ■ ^i))i whereas in- 
tersite attraction 14^ < introduce the PS state in which 
two NO phases coexists (with n^ = 2 and ?i_ = 0, i. e. 
electron droplets, PS[NO/NO]). For d = 1 the regions 
of the SS (or rjS, depending on the sign of /) and CO 
phases occurrence are extended, whereas a region of the 
NO (Mott) phase is reduced. The superconducting order 
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(in d = 1), similarly as in section III A is a short-range 



order whereas the CO is a long-range order (cf. also [46). 

The electron concentration n and chemical potential /i 
are (thermodynamically) conjugated variables in the bulk 
systems. However, one can fit the concentration rather 
than the potential in a controlled way experimentally. In 
such a case /i is a dependent internal parameter, which 
is determined by the temperature, the value of n, and 
other model parameters (cf. (pi)). Thus the obtained 
phase diagrams as a function of the concentration are 
quite important because in real systems n can vary in 
a large range. 

Although our model is (in several aspects) oversimpli- 
fied, it can be useful in qualitative analysis of experi- 
mental data for real narrow-band materials and it can be 
used for better understanding of various systems men- 
tioned below and in section IT] In particular, our results 
predict the existence of the phase separation (SS-NO) 
near the Mott state and describe its possible evolution 
and phase transitions with increasing temperature and 
a change of n (/l). 

Of course the PS instability is specific to the short- 
range nature of the interactions in the model. When 
the (unscreened) long-range Coulomb interactions are 
included the large-scale PS of charged particles is pre- 
vented and only a frustrated PS can occur (mesoscale, 
nanoscale) with the formation of various possible tex- 
tures |i7wg]. 

The electron phase separation involving SS is shown 
experimentally in several systems. In particular, for spe- 
cial cases of La2Cu04+5 and La2_2;Srj.Cu04_|_5, muon 
and superconducting quantum interference measure- 
ments suggest that the electron inhomogeneities move 
beyond local variations to form fully phase separated 
regions |S}|H]. Organic compounds also exhibit the 
superconductor-insulator phase separations as a result 
of the external pressure (e. g. quasi-one dimensional 
(TMTSF)2PF6 [lOj, (TMTSF)2Re04 Lllj) and fast cool- 
ing rate through the glass-like structure transition (e. g. 
K-(ET)2Cu[N(CO)2]Br fT^). Finally we mention the 
family of iron-pnictides, e. g. in Bai_a;Ka;Fe2-j,Se2, 
where mesoscopic PS between SS and insulating (mag- 
netic) phases have been observed [S^, T . 

The detailed analysis of RPA solutions will be explored 
further in subsequent publication (42]. It is of interest to 
analyze the effects of external magnetic field as well as 
the impact of density-density [46 , 50_ and magnetic [5T] 
interactions on the phase diagrams of model (fTl) . 
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Appendix A: Site-dependent self-consistent VA 
equations 

Within the VA the on-site interaction term is treated 
exactly and the intersite interactions are decoupled 
within the MFA (site-dependent): 

ptp] ^ {pt) Pj + (p7) pt - (pt) (pj) ■ (Ai) 

A variational Hamiltonian has the following form: 
^o = Y^ [Un,^n^l - nn, - 2x*pi - 2xipt + 

i 

+ x:A^ + X.A*], (A2) 

where Xi = ll]=ti ^ij^j' A* = (pf) and n^ = (n,). Hq 
can be diagonalized easily and a general expression for 
the grand potential J7 in the grand canonical ensemble in 
the VA is 

^ = -^ln{Tr[cxp(-/3i/o)]}, 

where /3 — l/{kBT). The average value of operator A is 
defined as 



Tr 



(^> = 



exp(-/3iJo)i 



Tr 



exp(-/3i?o) 



TrB means a trace of operator B calculated in the Fock 
space. 

The explicit formula for the grand potential obtained 
in the VA has the following form 



n 



^o. = 5]|x:a,+x.a: 



u 

2" 



/i- -ln(2Z,; 



(A3) 



where 



exp (/3C//2) + cosh (^^^/WT4\^\ 



and p — i-i — U/2. The expression for the average number 
of electrons at i-site is given by 



ni - 1 = (n,) - 1 = ^ (A4) 

and the average of the charge exchange operator at i-site 
is derived as 






(A5) 



so one has a set of 2N self-consistent equations to solve 



consisting of: (i) N equations in form of (A4| and (ii) 



A'' equations in form of ( |A5| ) (two equations for each site 
from N sites). 

The free energy is obtained as F = O + /x(iVe), where 

Ne = X^i ^i! (^e) — "-^ is ^^^ number of electrons in the 
system (n is defined by (Isl). 

The double occupancy Di of the site i is determined 
by the following equation: 



A = {ni^nii) = 



1 



1 exp(/3C//2) 






(A6) 



where n^ is defined by (A4| 



The solutions of the set ( A4 )-( A5 1 can correspond to 



a minimum, a maximum or a point of inflection of the 
grand potential (A3 1. In other words, the solutions of the 



set (A4|-(A5) fulfil the necessary condition for a mini- 
mum of il with respect to {Aijfl]^. To find the solutions 
corresponding to stable (or metastable) states of the sys- 
tem, one should find a minimum of fJ with respect to all 
A.. 

One can show easily that in the case of neglected spa- 
tial variations of the superconducting order parameter 
(i. e. A = Ai for every site i) Eqs. (A3l- (|A5| reduce to 
Eqs. ([7]) and (|9|-([l0]) obtained in section |llj 



Appendix B: The particle-hole transformation 

Let us introduce the particle-hole transformation 
J = J'^J^ on the alternate lattice defined as 



J: 



(-l)^'c,,, 



(for (T=t,;), 



where A^ = if i G A and \i — 1 Hi ^ B {A, B label two 
sublattices in the alternating lattice) . Then the operators 
appearing in Hamiltonian M are transformed to 



1 



ni^2- fii, 



(Bl) 



whereas the thermodynamic parameters and aver- 
ages: chemical potential, free energy, specific heat, 
superconducting-order parameter and double occupancy 
transform as follows 



Appendix C: Some analytical VA results 
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From the set (|9|-(ll0| one can obtain the relation be- 
tween the chemical potential and the electron concentra- 
tion in both homogeneous phases: 



fJ-NO 



(CI) 



In- 



1 + ^{n - 1)2 - n{n - 2) exp (-^U) 



^J'SS ^ loin- 1). 



(C2) 



The formula determining the second order SS-NO 
transition temperature Tss (the limit A — > of equa- 
tions (|9l-(10l) has the following form 



Pss 



U 



In 



sinh(^gg|/:igg| 
\n-l\ 



cosh {|3ss^J'Ss) 



(C3) 



where l3ss = '^/{ksTss)- 

Let us stress that Tss determined from the above equa- 
tion can be not only a transition temperature between 
stable phases but also a transition temperature between 
metastable phases (e. g. in the regions of PS states stabil- 
ity) or unstable phases. However, all second order tran- 
sition temperatures which were obtained in secti on |IV| by 
comparison of free energies fulfill ( C3 ) . From ( C3 1 one 
can obtain Tss as a function of n (or p, - cf. (C2|). In 
the limiting cases it yields 



keTssiU ^ -(») = 2kBTss{U = 0) = 



2/o(n- l)ln"^ 



(C4) 



The TCP-line is described by the following equation: 

-1 



tanh(/3T/o|n- 1|) 



1 



1\[3tIo 



1 



(C5) 



Projections of TCP points on the kBT/lQ-n/Io plane are 
shown in figure |2] by dashed-dotted curve. 



/i(2 -n) = -n{n) + U, ( fl{2 - n) = 


--Kn) ) (B2) 


fi2-n) = fin) + Uil-n), 


(B3) 


c{2 — n) = c{n), 


(B4) 


A(2-n) = -A*(7i), 


(B5) 


D{2 -n) = D{n) + 1 - ?i. 


(B6) 



respectively. From any characteristics of the considered 
model in the range < n < 1 we can find the behaviour 
of the system in the range 1 < n < 2. In particular, the 
phase diagrams are symmetric with respect to half-filling. 
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